NACA TN 3280 Z L S 8 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

TECHNICAL NOTE 3280 


a 


tr 

tr 

ru 

JT 

UJ 




ELECTRICAL ANALOGIES FOR STIFFENED SHELLS 
WITH FLEXIBLE RINGS 
By R. H. MacNeal 
California Institute of Technology 



Washington 
December 1954 


^ n r. 


^ ' P. \ 


lU^Lnrrrr 


c *v. 

L 


TECH LIBRARY KAFB, NM 


TECH LIBRARY KAFB, NM 


NATIONAL ADVISORY COMMITTEE FOR AERONAUTIC 


0bb24: 


11 


TECHNICAL NOTE 3280 

ELE CTRICAL ANALOGIES FOR STIFFENED SHELLS 
WITH FLEX IB LE RINGS 


By R. H. MacNeal 


SUMMARY 


Structural theory and analogous electrical circuits are developed 
for stiffened shells with flexible rings. By assumption, the forces 
that a shell (consisting of stringers and skin) and a ring can exert on 
each other sure directed along their common line of intersection. As a 
consequence the shell and the rings can be treated separately. 

First an electrical analogy is developed for a circular shell with 
a straight axis and variable radius . This analogy is extended to non- 
circular cylinders . Next an electrical analogy is derived for rings with 
variable radii of curvature; a simplified circuit for circular rings is 
also presented. The simplifications that occur when the rings are assumed 
to be rigid are discussed. Finally results are given for two sample 
problems solved on an analog computer . The second problem concerns a 
cantilever conical shell and illustrates the manner in which the shell 
and ring circuits are interconnected. 


INTRODUCTION 


The type of shell considered in this paper has an elongated shape 
and consists of a thin skin supported by stringers and rings. In 
analyzing such shells it is the nearly universal practice to replace the 
elastic supporting rings by rigid bulkheads in order to simplify the 
analysis. This will not be done in this paper. 

The means of analysis to be used in this paper is an electric analog 
computer of the "direct analogy" type. Any complicated system, if it is 
to be analyzed on such a computer, must have its equations formulated in 
a very special way. Essentially one seeks for laws of equivalence between 
the system being analyzed and a lumped constant electrical network. The 
basic laws of equivalence between the equations of elasticity .and the 
equations of an electrical circuit are well known. In fact there are two 
alternative sets of laws depending on whether force is made analogous to 
current or to voltage. If the former alternative is chosen the laws of 
equivalence are: Force is analogous to current, displacement is a n alogous 
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to voltage, Hooke's law is analogous to Ohm's law, equations of equilib- 
rium are analogous to Kirchboff 's law for the sum of currents entering a 
node, and the equations concerning the compatibility of strains are 
analogous to Kirchboff 's law for the voltages around a loop. 

However much comfort these basic laws of equivalence may give they 
are us uall y insufficient to determine the form of a lumped- constant 
electrical network that is analogous to a given structure. For one thing 
elastic structures are continuous rather than "lumped," and some means 
must be found for replacing the given continuous elastic structure by an 
Idealized lumped one before an electrical analogy can be found. This 
"lumping" consists either of replacing the differential equations 
governing the structure by finite-difference equations, or of employing 
other devices such as concentrating normal-stress -carrying area into 
equivalent flanges and shear-carrying area into equivalent panels. In 
the analysis of stiffened structures this latter approach is reinforced 
by the fact that much of the structure is in fact so concentrated. In 
this paper both of the methods mentioned will be used. 

In deriving an electrical analogy for an elastic structure an effort 
should be made to preserve a one-to-one correspondence between the prop- 
erties of the electrical circuit and the properties of the idealized 
structure. This correspondence means , for example, that the current in 
resistor A is equal to the force in flange A' multiplied by a scale factor, 
or that the voltage at node B is equal to the vertical displacement at 
panel point B’ multiplied by a scale factor. If such correspondences Eire 
preserved, the analog computer can be made a useful tool for designing 
as well as for analyzing structures . If a change in the cross-sectional 
area of a single flange corresponds to changing the value of a single 
resistor and if currents can be easily and directly converted into 
internal forces, then design changes can be made very rapidly and their 
effects instantly determined while the problem is set up on the analog 
computer. In the present paper these correspondences sire rigidly 
preserved. 

Another advantage of the close one-to-one correspondence of the 
electrical analogy and the idealized structure is that it enables the 
structural engineer, who is usually uninstructed in electric-circuit 
theory, to understand the operation of the analog computer and to use 
it himself after a period of indoctrination. It has even been found 
that structural engineers may be aided in their understanding of structures 
by using some of the concepts of electric-circuit theory. This naturally 
applies, a fortiori, to the electrical engineer. 

The present paper is intended as a step in the development of an 
analog-computer method that is generally applicable to the solution of 
aircraft structural problems. At present, analogies exist in the tech- 
nical literature for beams, frameworks, flat sheets, the bending of 
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plates , and the tending of platelike multicell shells. Structures 
combining components of the above types can be analyzed by combining 
their electrical analogies . Consequently it is at present possible to 
analyze, a great many practical aircraft structures. 

Some of the previous papers that have a direct bearing on the sub- 
ject of the present paper should be mentioned. 

In 19^4 Kron published a paper containing electrical analogies for 
the general three-dimensional elastic-field problem and, as subcases, 
analogies for the plane-stress and plane-strain problems (ref. l) . In 
a companion paper. Carter worked the plane-stress problem for a deep 
cantilever beam (ref. 2). 

More recently analogies have been developed for thin multicell 
shells having a horizontal plane of symmetry by using an equivalent plate 
theory (ref. 3 ) • 

In 1951 Goran published a paper containing an electrical analogy 
for stiffened elastic shells (ref. !<-) . The shells were assumed to be 
conical and to be supported by rigid bulkheads. Although the stringers 
were not assumed to be parallel, the panels were assumed to be nearly 
rectangular in shape, Goran used a minimum energy principle in deriving 
the equations from which he developed the electrical analogy, in contrast 
with the method of difference equations used in this paper. 

The present investigation was conducted at the California Institute 
of Technology and has been made available to the national Advisory 
Committee for Aeronautics for publication because of Its general interest. 


SYMBOLS 


■A-s 

E 


Fn 

f b ,f s ' 

Fst> Fg-t’ 


cross-sectional area of stringer 
Young's modulus 

external load normal to ring per radian of <f> 
external load tangential to ring (or skin) per radian of 0 
shear force in ring 
forces in stringer 

tangential forces in panel parallel to t 
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As,At,A0 


axial force In ring 

tangential forces in panel parallel to s 
shear modulus 
thickness of skin 

moment of inertia of ring cross section 
tending moment in ring 
applied vertical load 

radius of circular shell; ' radius of curvature of ring; 
distance to point in rigid ring 

coordinate parallel to stringer 

coordinate perpendicular to stringer and parallel to ring 

displacement normal to axis of ring 

displacement parallel to s 

displacement parallel to t 

coordinate related to ut hy transformation 

horizontal displacement of rigid "bulkhead 

horizontal direction 

vertical displacement of rigid bulkhead 
vertical direction 

angle between two adjacent stringers; angle between tangential 
displacement and line to a point in a rigid bulkhead 

shear strain of panel 

angle by which direction of s-axis is changed because of 
translation in t-direction 

difference operators in s-, t-, and ^-directions 

increments in a, t, and 0 
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A0' = 2 sin 

j 1 5 angle between horizontal line and normal to center line of 

undisplaced ring 

$ angle between stringer and axis of shell 

8 rotation about axis of shell, usually in ring 

DERIVATION OF AN ANALOGY FOR A CIRCULAR NONCYLIKDRICAL 

STIFFENED SHELL 


A sketch of a circular noncylindrical shell is shown in figure l(a) . 
This shell consists of circular elastic rings to which stringers and a 
thin skin covering are attached. The rings are spaced a finite distance 
apart in planes perpendicular to the axis of the shell, which is assumed 
to be straight. The stringers lie in planes perpendicular to the rings 
and are assumed to carry axial forces only. The skin is divided into 
panels by the intersection of the stringers and rings; these panels are 
assumed to carry shearing forces only, the effective normal-stress - 
carrying area of the sk in having been included in the cross-sectional 
areas of the stringers and rings. The radius of the shell may vary in 
any manner along the axis of the shell, but the radius of curvature of 
the stringers is assumed to be large compared with the radius of the 
shell. The angle between the axis of the stringers and the axis of the 
shell need not be small provided that the rings and stringers are spaced 
close enough together so that the shear panels are approximately 
rectangular in shape. 

The stringers and skin will be treated separately from the rings. 
For the analysis of the stringers and skin it will be shown that the 
rings can be represented by tangential external forces applied to the 
skin along the lines of intersection of the skin and the rings. 

The equilibrium and force -displacement equations will be derived 
for the skin and stringers and an electrical circuit satisfying these 
equations will be constructed. Then in the next section the equations 
for an elastic ring will be written and a circuit satisfying these equa- 
tions will be constructed. As a final step the two circuits will be 
connected together to give the electrical analogy for the whole shell. 
Orthogonal coordinates in the surface of the shell parallel and perpen- 
dicular to the stringers as shown in figure l(a) will be employed in the 
analysis of the skin and stringers. An enlarged portion of the shell is 
shown in figure l(b) . This figure shows the points where displacements 
parallel to the axis of the stringers Ug and displacements parallel to 
the axis of the rings u-^ are defined. 
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Figure 2 sliovs a portion of the skin "between two adjacent rings 

with its midpoint on a stringer. The total axial forces carried "by the 

stringer at points where the stringer passes over two adjacent rings 

are and F 0 ' . The total tangential forces acting in the s -direction 

s s 

on sections passing through the centers of two adjacent shear panels are 
F ts and F^’. 

The force in the stringer is continuous at the point where it passes 
over a ring because of the following assumptions: 

(a) Thp ring can exert only forces which lie in its own plane 

(b) The stringer cannot support bending loads 

(c) THp change in direction of the stringer at the point where it 
passes over a ring is negligibly small compared with the curvature of 
the ring. J- 

Consequently the ring can exert forces only in the t -direction, and these 
forces can be treated as applied forces in the analysis of the skin and 
stringers . 

In figure 2, a/2 is the angle between the stringer and the line 

of action of the shear force F-fcs ' * This angle is approximately equal to 

(l/2)ijr where i}r is the angle between the axis of the shell and the 
axis of the stringer and A0 is equal to the angle subtended by a segment 
of ring between two adjacent stringers. The cosine of the product of 
these angles will be assumed to be equal to 1. Hence the equilibrium 
equation for forces in the s-direction is: 


V 


- F s + F ts’ 


- F ts = 0 


( 1 ) 


Using difference-equation notation, equation (l) can be written 


Vs + Vts = 0 


( 2 ) 


This notation will be used in the remainder of this paper. 


It is not difficult to relax this assumption to permit an appreciable 
change in direction of the stringer and a consequent radial load on the 
ring. The resulting electrical circuit contains an additional transformer 
at each point of intersection of stringer and ring. 
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Figure 3 shows a portion of the skin between two stringers with its 
midpoint on a ring. The total tangential forces acting in the t-direction 
on sections pass in g through the centers of two adjacent shear panels are 
F a t find Fg-fc ' . The tangential force exerted by the ring on the skin per 
radian of 0 is ffc. The lines of action of the forces F^ s and -^ts 

intersect the axis of the shell, so that the equation of equilibrium for 
moments about this axis is 


As (rFg-t) - rA0ft = 0 


(5) 


The total forces acting on perpendicular planes passing through the center 
of a shear panel are F-^g and F s t. From equation ( 3 ) it is seen that 

the shear stress cannot be uniform in the s-direetion across the surface 
of a panel if r is not constant. Hence the assumption that the panel 
carries only shear stresses is incorrect. The secondary normal stresses 
required by equation ( 3 ) are ignored. 

It will be assumed that the variation of shear stress is linear 
across the surface of the panex, so that the value at the center of the 
panel is equal to the average along a line in either the s- or t-direction. 
Then a relationship between Fts and F s t may be obtained from the 

equilibrium equation of a small element at the center of the panel (see 

fig. 4): 


1 8t - I 68 <« 


Fts - Fst gj (5) 

The equilibrium of the portion of the shell shown in figure 3 for 
forces in a direction parallel to t at the center of the section may be 
demonstrated, if desired, by means of equations ( 3 ) and ( 5 ) . 

The force-displacement equation for the stringers is quite easily 
written. The axial d isp lacement of the stringer Ug is defined at the 
midpoints between adjacent rings as shown in figure l(b) while the axial 
force (positive for tension) is defined at points where the stringer 
passes over the rings as shown in figure 2. Ass uming the variation in axial 
force to be linear between adjacent points where Ug is defined, from 

Hooke’s law 
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F 6 “ ^ (Ab%) 

where Ag is the cross-sectional area of the stringer. 

The relationship between shear force and displacements for the shear 
panels is less easily written because the relationship between shear 
strain and the displacements - is complicated. The displacements in the 
s- and t-directions are defined at the midpoints of the sides of the 
shear panel as shown in figure l(b) . The shear strain of the panel is 
defined as the distortion of the angle between two lines passing through 
the midpoints of the sides . In computing this angle care must be taken 
to eliminate apparent distortion due to rigid body rotation about the 
axis of the shell. From figure 5* the shear strain 


7 = 


y l + 


(T) 


where 


, = fVs 

1 At 


( 8 ) 


-v - ^ u t 
7 2 . ~ As 


6 


(9) 


and 8 is the angle by which the direction of the s-axis, drawn through 
the center of the panel, has been changed because of translation parallel 
to the t-aodLs. For small displacements this angle is 


( 10 ) 


where ljr is the angle between the axis of the shell and the s-axis. 
From figure 2 It can be seen that 


8 = — sin tjr 
r 


= sin ijr 
8s 


(H) 
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Combining equations (9), (lO), and (ll). 



Vt _ H dr 
As r 8s 


( 12 ) 


Assume that the first term on the right can he replaced h y the 
equivalent partial derivative evaluated at the center of the panel 


_ ^t ^ 8r r _8_K\ 

8 s r 8 s 8 s \r j 


(13) 


Replace this equation by its finite difference equivalent and obtain 
from equations ( 7 ), ( 8 ), and ( 13 ): 


7 = 


At + 


(HO 


The shearing strain is related to the total tangential force acting 
on the panel by the following equation 


7 = 


_L 

Gh As 


(15) 


where h is the thickness of the skin. Hence the force-displacement 
equation for the panel is 


'*■ - 011 § 




( 16 ) 


The equations that are essential for the construction of an electrical 
analogy are summarized below. The equilibrium equations are: 


Vs + ^t-^ts — ® 


(17a) 


A B (rT st ) - Atf t = 0 


(17b) 
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The force-displacement equations are: 


F 


s _ 



F ts = Gh 
158 At 


At(us) + ^ 


(174) 


(I7e) 


These equations have heen derived hy assuming the normal-stress- 
carrying area of the shell to he concentrated in stringers. They can 
also he derived from the differential equations for a membrane shell of 
revolution (ref. 5 ) "by replacing differential operators hy finite dif- 
ference operators. This is an important fact because it extends the 
applicability of the equations to certain unstiffened shells. 

In the electrical analogy forces are analogous to currents and 
displacements are analogous to voltages. The complete circuit is shown 
in figure 6. This circuit consists of two separate parts. In one part 
the voltages to ground are the displacements u g while in the other 

part the voltages to ground are the rotations u^/r. The two circuits 

are coupled together hy means of ideal transformers . Transformer coils 
which are coupled together are indicated hy circled numbers. Points at 
which each one of the above equations are satisfied are indicated hy 
letters in the circuit . Equation (l7&) is satisfied hy the currents 
entering a node of the Ug circuit. Equation ( 17 b) is satisfied hy the 
currents entering a node of the u-fc circuit. Equation ( 17 c) is satisfied 

hy the currents flowing in the windings of a transformer whose turns 
ratio is rAt/As . Equation (l7d) is satisfied hy a resistor whose value 
in ohms is As/EA S . In equation (l7e) the increment in Ug in the 
t-direction is added to a fraction of the increment in u-^/r in the 

s-direction. This addition is accomplished hy the same transformer which 
satisfies equation ( 17 c) . The sum of these terms is the voltage across 
a resistor whose value is At/GhAs and through which a current equal 
to F.fc S flows . 

Two observations can he made concerni n g equations ( 17 ) and the 
resulting circuit. If r does not depend on s , the equations are those 
of a cylindrical shell, and r may he removed from Inside the difference 
operators. Hence the equations of a noncylindrical shell have the same 
form as the equations of a cylindrical shell if rotation about the axis 
u^/r (rather than tangential displacement) and torque about the axis 

rF 8 -fc (rather than tangential force) are used as variables. The rotation 

and torque about the axis are the natural variables to use in deriving 
the equations of a circular noncylindrical shell. 
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Equations (17) also apply with slight modification to the skin and 
stringers of a noncircular cyl indrical shell. In this case r depends 
on t rather than on s, hut this dependency on t will not enter 
into the derivation of the equations for the skin and stringers. Hence 
for a noncircular cyl indrical shell, r can he eliminated from equa- 
tion (l7e) and can he divided out of equation (17b) to give 


Ag(Fst) " A^t = 0 


( 18 ) 


Since ft is the external tangential force per radian of 0 , I s 

the total external tangential force per hay in the t-direction. 

In figure 6 parts of the ut circuit corresponding to different 

values of t are not shown connected. The interconnection is accomplished 
hy means of the currents Atft which are the reactions of the rings. 

The voltages at the points where these currents are inserted are con- 
strained to he equal to the corresponding values of u^/r for the rings. 

A complete circuit for a shell, including the circuits for elastic rings, 
will he shown later. If the rings are assumed to he rigid, the circuits 
for the rings "become quite simple as will he shown. 


DERIVATION OF AN ANALOGY" FOR AN ELASTIC RING 


An electrical analogy for the bending of a ring in its own plane 
will he derived in this section. In the discussion of shells supported 
hy rings it is usually assumed that the rings are perfectly free to warp 
out of their own planes . An electrical analogy for a circular ring 
deforming perpendicular to its own plane has been derived hy Russell 
(ref. 6), hut this effect will not he considered here. 

It will also he assumed that the effects of shearing stiffness and 
axial stiffness of the ring are small so that these effects can he 
ignored. This assumption is made in order to simplify the discussion; 
these effects can he included in the electrical analogy if desired. In 
addition the eccentricity of the ring, that is, the distance between the 
neutral axis of the ring and the point of attachment to the skin, will he 
assumed to he zero. This important effect can also he included in the 
electrical analogy. 

An element of ring is shown in figure 7 * The displacement quantities 
to he used in the analysis of the ring are the normal and tangential 
displacement of the center line and the rotation of the normal to the 
center line. The independent position variable 0 is the angle between 
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a horizontal line and the normal to the center line of the unloaded ring. 
The external tangential and normal loads per radian of <f> applied along 
the center lin e of the ring are f-fc and f n . Distributed moment loads 

will not be considered. It will be noted in figure 7 that the radius of 
curvature is not assumed to be constant. 

The loads, internal forces, and displacements of the above described 


ring satisfy the following 
equilibrium equations are: 

six first-order differential equations. 

The 

• 

dFn 

— a = F+ - f - 
d0 t n 

(19) 


dF + 

(20) 

• 


(21) 

The stress-strain and strain-displacement equations are: 



d0 _ Mr 
d0 El 

(22) 


u 

CD 

+ 

(23) 


■ejl 

11 

£ 

(24) 


These equations will be replaced by the corresponding first-order 
difference equations. In so doing central difference equations will be 
•used so that the quantities appearing be hind the derivative symbols in 
the above equations are defined at values of midway between the values 
of $ at which the undifferentiated quantities are defined. For example, 
equation ( 19 ) may be approximated by the following equation: 
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The same difference-equation notation can be used for this equation 
as was used in the previous section: 

= (F t - fn) A 0 ( 26 ) 

Position subscripts are not required in this equation. The quantities 
u n , Ft, ant ^ M are defined at the same points and these points are mid- 
way between the points where 0 , u^, and F n are defined. 

Equations (20) to ( 2k ) could be replaced by simple difference equa- 
tions in the same manner that equation (19) was replaced by equation (26). 
However, it can be demonstrated that the following difference equations 
are more accurate. They give exactly correct results for a segment of 
ring which is rigid, has constant radius of curvature, and is uniformly 
loaded. In other words, for any complete ring with constant radius of 
curvature , they give correct results for rigid body displacements and 
the statically determinant parts of the internal forces. These statements 
require lengthy proofs and instead of the proofs being given they will be 
partially demonstrated later by means of an example. 

The equilibrium equations are: 

= (*t - f n) (2 Sto (27) 

y t - -( r n + f t) (2 ^ (28) 



The stress-strain and strain^displaceinent equations are: 

V - ^ M (50) 

= H + r 0 ) ( 2 sin if 


(31) 
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= _u n^ sln ^ (52) 

For small values of £0, 2 sin (£0/2) is approximately equal to £0 

and the above equations approach the corresponding simple central- 

difference equations. The chord subtending the arc £0 is r^2 sin 

An electrical circuit which identically satisfies the above equations 
is shown in figure 8. In this circuit displacement quantities axe voltages 
to ground and the loads and internal forces are currents. Equation ( 29 ) 
is satisfied by the currents entering a junction in the upper circuit of 
figure 8. Equations ( 27 ) and ( 28 ) are similarly satisfied by the currents 
in the middle and lower circuits of figure 8. Ideal transformers axe used 
to produce currents flowing into the junctions of one circuit that axe 
proportional to the currents flowing in the "main line" of another circuit. 

Equation ( 30 ) expresses Ohm's law for the drop in voltage between 
successive nodes of the upper circuit. Equation ( 31 ) is satisfied in the 
main line of the middle circuit by means of ideal transformer coils which 
insert voltages in the line proportional to the voltages to ground in the 
other two circuits. Equation ( 32 ) is similarly satisfied in the lower 
circuit. Each transformer is instrumental in the satisfaction of two 
equations, an equilibrium equation and a strain-displacement equation. 

The circuit of figure 8 employs three transformers per cell. If the 
radius of curvature of the ring is constant (fig. 9(a)), or if the ring 
can be divided into segments containing several cells for each of which 
the radius of curvature is constant, a circuit requiring only two trans- 
formers per cell can be used. Since r is now assumed to be constant. 


equations ( 27 ) to ( 32 ) can be rewritten as follows (where 
has been abbreviated by £0 ' ) : 

2 sin £0/2 

Afn = (*t r )^r- - f n^’ 

(27a) 

^( F t r ) = -F n (r£0') + ( f t r )^' 

( 28 a) 

Af. = -F n (r£0') 

(29a) 




( 50 a) 



MCA TN 3280 


15 


- (-9)j rA^ 1 (31a) 

- u n(^) ^2a) 

In this form of the equations uj./r and F^r replace u^. and F-f. 

as variables as in the case of the skin and stringers for a noncylindrical 
shell. The purpose of this manipulation has been to get equation (31a) 
into the form shown, where the increment in u n is proportional to the 

difference of the other two displacement quantities. This equation can 
he satisfied by a single transformer whereas equation (31) required two 
transformers. A circuit satisfying equations (27a) to (32a) is shown 
in figure 9(h) . In this figure the transformer coils connecting the 
upper and middle circuits are coupled to the transformer coils in the 
main line of the lower circuit. T hi s remote coupling is indicated by 
circled numbers. The currents corresponding to the loads are not shown 
in this circuit. 


CIRCUITS FOR RIGID RINGS 


It Is frequently possible to assume that some or all of the rings 
supporting a shell are rigid in their own planes without serious error 
in the analysis of the shell. This assumption greatly simplifies analyt- 
ical solutions of shell problems and sometimes eliminates a great deal 
of the equipment required in an analog-computer solution. Since it Is 
assumed here that only the tangential displacement of the ring is impor- 
tant in the analysis of shells, this is the only coordinate that need be 
represented at points around the periphery of a rigid ring. 

The position of a rigid bulkhead is determined by the displacement 
of one of its points in two perpendicular directions and by the rotation 
about an axis perpendicular to Its plane (as shown in fig. 10(a)). The 
tangential displacement at points on the periphery can be computed from 
these three quantities 


u-fc = Y sin <f) + Z cos ft - 9r sin a 


(33) 


In an analytical solution Y, Z, and 0 are unknown quantities. They 
are usually regarded as Lagrangian multipliers and equation (33) i 0 
regarded as an equation of constraint. In an electrical analogy this 
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equation of constraint can "be satisfied "by a network of transformers, a 
general form of which is shown in figure 10(b) . This network also 
satisfies the equilibrium equations of the rigid ring. Applied loads 
in the y- and z-directions and applied torque are inserted as currents 
into the network as shown. Since, in general, this circuit requires 
three transformers for each tangential displacement it has no advantage 
over an elastic ring circuit in which the resistors corresponding to the 
bending stiffness of the ring are set equal to zero. However, in prac- 
tice, one or more of the terms in equation (33) may t>e equal to zero 
because either the unknown or the coefficient may vanish. 1 In such cases 
the transformer network may be quite simple. 

For example, in a shell with a vertical plane of symmetry, loaded 
symmetrically with respect to this plane. 


u-fc = Z cos 0 


(34) 


This equation requires one transformer for each value of u^. It is 
possible to introduce further simplification by replacing ut by a 
variable which depends on 0. Let 


u t = cos 0 


(35) 


then 


u t = Z 


(36) 


If this change of variable is now introduced into the equations of the 
skin and stringers (eq. (17)), it will be found that the form of the 
equations will remain unchanged and that the only effect will be a change 
in the turns ratio of the transformers coupling the Ug and uj. circuits. 

Equation (36) is then satisfied by connecting together all the tangential 
displacement nodes in any one ring. 

As another example consider a rigid ring (or bulkhead) which is very 
thin in one direction (see fig. 11(a)). In this case it may be assumed 
that the angle between the y-axis and the top and bottom surfaces of the 
shell is small everywhere except at the ends, where there are closing 
vertical segments. Furthermore it is usual to assume that Y, the 
displacement para! lei to the long direction, is zero. In this case 
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u-fc = -0r sin a 
ut = Z - (0c/2) 
^ = -Z - (0c/2) 


top and bottom surfaces 
at left end 
at rlglrfc end 


•\ 




j 


( 37 ) 


Here again scale factors can be introduced to simplify the equations: 


Ut = - 




= 0 


r sin a 
Ufci = u ti = Z - (0c/2) 


top and bottom surfaces 


at left end 


> ( 38 ) 


= -U-j. = Z + (0c/2) at right 


end 


These equations are satisfied by the simple network of figure ll(b) . 

With. this circuit for a rigid bulkhead and the circuit for the 
stringers and spars (fig. 6), two spar box wings with unsymmetrical top 
or bottom surfaces can be analyzed. The extension to multispar wings is 
simple and direct. 


SOLUTION OF PROBLEMS 


The Cal-Tech analog computer (ref. 7 ) was used for the solution of 
two problems in connection with the preparation of this paper. This 
computer consists essentially of a storehouse of electrical parts which 
contains , among other things, the resistors and high-quality transformers 
required in the solution of stress-analysis problems. 

The first problem was the analysis of the simple circular ring 
shown in figure 9 (a.) subjected to two opposing concentrated radial loads. 
Because of the symmetry of the ring and its loads a quadrant of the ring 
can be substituted for the whole if proper boundary conditions are applied 
-at the ends of the quadrant. At 0 = 0 the proper boundary conditions 
are that 0, ut, and F n equal zero. At 0 = rt/2 the proper boundary 

conditions are that 0 and ut equal zero and that F n equals P, one- 
half of the applied load. 
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In figure 9(1), tie quadrant of ring has "been represented, by a 
circuit conta ining four cells . The boundary conditions have been 
satisfied by setting the corresponding electrical quantities equal to 
zero at the two ends. 

The results of the analysis are presented in table I in dimensionless 
form. An exact analysis of the problem using differential equations is 
compared with the analog-computer solution. In addition an exact solution 
of the difference equations governing the electric circuit is shown. It 
will be seen that the differential-equation solution and the exact solu- 
of the difference equations give identical results for the internal shear 
and inter nal axial force. This was to be expected since these quantities 
are static ally dete rmina te in the problem investigated. The other quan- 
tities show errors due to finite-difference approximation. A comparison 
of the difference-equation solution and the analog-computer solution 
shows errors in the computer solution of the order of 1 or 2 percent, 
which is fairly typical of the results customarily obtained with the 
Cal-Tech analog computer . 

The second problem was the analysis of a conical shell supported by 
circular elastic rings. As such it provides an example illustrating the 
manner in which the shell and ring circuits are interconnected. The 
structure, which has 3 elastic rings and 14 stringers, is shown in fig- 
ure 12(a), and specifications- for the structure are given in table II. 

This structure is supposed to resemble the aft portion of an aircraft 
fuselage. The structural weight is divided approximately equally between 
the skin, stringers, and rings. The rings are somewhat stiff er than those 
employed customarily in fuselage construction. The number of stringers 
nnd the number of rings in the structure are much fewer than the number 
that would be employed in an aircraft fuselage, so that each stringer of 
the structure represents several stringers in the fuselage, and the stiff- 
ness of intermediate fuselage rings is included in the stiffness of the 
three main rings shown. 

The shell is subjected to symmetrical concentrated vertical loads 
applied to the ring at the small end of the shell and these loads are 
reacted at the large end, which is built into a rigid wall. Because of 
the symmetry of the structure and the applied loads, only a quarter of 
the shell need be considered if appropriate boundary conditions are applied. 
This part of the shell has been given a two-coordinate numbering system 
for the identification of points in the structure . For example, point 42 
refers to a point at the intersection of ring number 4 and stringer number 2. 

The electrical analogy for the structure is also shown in figure 12. 

This circuit consists of three parts, the u s and u^ circuits and the 

ring circuits. In figure 12(b) only one ring circuit is shown since the 
other two have an identical appearance. The connections between the ring 
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circuit shown and the u-fc circuit are indicated by the circled letters a, 
b, and c. Currents corresponding to the interaction forces between the 
skin and the rings flow through these connections. In figure 12(a) a 
cutout is Indicated in the middle bay. Electrical parts corresponding 
to the cutout (shown dotted in figs. 12(c) and 12(d)) are removed when 
the cut panel is removed. 

The boundary conditions at the vertical plane of symmetry are that 
the shear stresses in shear panels intersected by this plane are zero 
and that in the ring circuits 0, u-fc, and F n are al 1 zero except when 

vertical loads are applied to the ring in the plane of symmetry, in which 
case F n is equal to P, one-half of the applied load. The boundary 
conditions at the horizontal plane of symmetry are that the displacement 
in the s-direction Ug is zero and that in the ring circuits u^, M, 

and F-fc are zero except when vertical loads are applied to the ring in 
the horizontal plane of symmetry, in which case F^ is equal to -P. 

The boundary conditions at the large end of the shell are that both Ug 
and U£ are equal to zero. At the small end the boundary conditions are 
conveniently expressed as the absence of applied loads except as indicated. 
All of the above boundary conditions have been satisfied in the electrical 
circuit by means of short and open circuits. Short circuits are used to 
set displacements equal to zero while open circuits are used to set internal 
forces equal to zero. 

The electrical circuits shown in figure 12 use 27 resistors and 
27 transformers. In calculating the values of the transformer turns 
ratios and resistors to be used in the circuit it is necessary to make 
use of scale factors. This aspect of the problem has been omitted in 
the present discussion in order that the electrical quantities shown in 
the circuit diagram may have a direct significance in terms of mechanical 
quantities. A brief discussion of the scale-factor methods employed with 
the Cal-Tech analog computer is given in the appendix of reference 8. 

The structural and loading conditions that were investigated in this 
problem are given in table H. It was desired to investigate the effect 
of the following things on the distribution of internal forces: The 
effect of the stiffness of the rings, the effect of the location of the 
vertical load on the end rings, and the effect of a cutout in the middle 
bay. 


The results of these investigations are given in table II. The 
quantities recorded according to the numbering system previously dis- 
cussed are the stringer forces F s , the panel shears F^ s , the ring 

bending moments M, and the vertical displacements. The vertical dis- 
placements have been obtained by vectorially combining the tangential 
and normal displacements of the rings . The tabulated results are subject 
to experimental errors and when four significant figures are given the 
fourth figure is entirely unreliable. 
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The main conclusion to "be drawn from these results is that for the 
shell tested and the loads which were applied the stiffness of the rings 
has a rather small effect on the distribution of internal forces. The 
only change made between conditions (l) and ( 2 ) was that in condition ( 2 ) 
the rings were made five times as stiff as in condition (l) . It will be 
seen that the stiffness of the rings has very little effect on the 
internal forces except between the second and third rings. However, the 
distortion of the third ring in condition (l), as given by the vertical 
displacements, is significant. In another condition, the results of 
which are not tabulated, the first and second rings were stiffened by a 
factor of 10 while the third ring retained its normal stiffness. The 
di f ference in the results between this condition and condition (l) was 
negligible . 

The results of condition ( 3 ) indicate that the effect of ring stiff- 
ness is considerably less when the applied vertical loads are in a 
horizontal plane than when they are in a vertical plane. The results of 
conditions (4) and ( 5 ) indicate that even moderately small symmetrical 
cutouts produce a severe redistribution of the internal forces and that, 
in this case, the effect of ring stiffness is important if accurate 
results axe desired. 


California Institute of Technology, 
Pasadena, Calif., November 25, 


1953 
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TABLE I.- SOLUTION OF CIRCULAR-EING PROBLEM WITH 
FOUR FINITE-DIFFERENCE CELLS PER QUADRANT 


Quantity 


Differential 

Difference 

Analog 

equations 

equations 

computer 

j* 

22.5 

0.3827 

O.3827 

0.387 

X n 

45-0 

.7071 

.7071 

.702 

P 

67.5 

.9239 

.9239 

.922 


11.25 

0.9808 

0.9808 

0.982 

F t 

33-75 

.8315 

.8315 

.824 

p 

56.25 

.5556 

.5556 

.566 


78.75 

.1951 

.1951 

.210 


H.25 

0.3442 

0.3401 

0.340 

M 

33-75 

.1948 

.1907 

. .192 

Pr 

56.25 

-.0810 

-.0851 

-.083 


78.75 

-.4415 

-.4456 

-.447 

_SL 

22.5 

0.1327 

0.1335 

0.1338 

Pr 2 

45-0 

.2071 

.2084 

.2064 

El 

67.5 

• 1739 

.1750 

.1703 

^t 

22.5 

0.0487 

0.0514 

0.0515 

Pr 5 

45-0 

.0706 

.0747 

.0756 

El 

67.5 

.0515 

.0548 

.0559 


11.25 

-0.1271 

-0.1318 

-0.1310 

Pr 3 

33-75 

-.0573 

-.0596 

-.0586 

56.25 

.0493 

.0509 

.0535 

El 

78.75 

.1349 

.1405 

.1435 
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TABLE H.- RESULTS FROM STUDY OF CONICAL SHELL 


Specifications for shell: Number of stringers, 14; cross-section area 

of stringers, 0.60 sq. in.; thickness of skin, 0.050 in.; Young's 
modulus, 10.4 x 10^ psi; and shear modulus, 4.0 x 10^ psi. P = 1 kipT] 


Panel 

Condition 

(1) 

(2) 

(3) 

(4) 

(5) 

point 


(a) 

(b) 

(c) 

(a) 

(e) 



Stringer forces, F s , 

kips 


2 

-0.945 

-0.945 

-0.942 

-0.330 

-0.577 

4 

-1.708 

-I.708 

-1.708 

-1.550 

-1.632 

6 

-2.118 

-2.050 

-2.124 

-2.475 

-2.300 

22 

-.688 

-.711 

-.722 

.767 

.688 

2 k 

-1.282 

-1.282 

-1.291 

-1.409 

-1.321 

26 

-1.612 

-1.594 

-I.587 

-2.124 

-2.180 

42 

-.308 

-.380 

-.452 

-1-7® 

-1.789 

44 

-.676 

-.714 

-.745 

-.536 

-.650 

4 6 

-1.002 

-935 

-.882 

-.445 

-.326 


Panel shears, F-(. a , kips 

11 

1.206 

1.195 

1.191 

1.593 

1.726 

13 

• 954 

.956 

.963 

.509 

.469 

15 

.513 

• 531 

.536 

.372 

.155 

31 

1.596 

1-550 

■ 1.522 

0 

0 

33 

1.217 

1.231 

1.254 

2.537 

2.457 

35 

.618 

.661 

.707 

1.687 

1.838 

51 

1.985 

2.046 

2.080 

2.740 

2.789 

53 

1.685 

1.662 

1.632 

•994 

1.032 

55 

1.013 

.937 

.880 

.449 

-334 



Ring bending moments, M, 

in-kips 


22 

-0.07 

-0.05 

0.02 

3-58 

4.08 

24 

-.03 

-.02 

.01 

-96 

1.28 

26 

.06 

.02 

-.02 

-2.40 

-2.84 

42 

.20 

.05 

-.08 

-3.72 

-3.88 

44 

.10 

.01 

-.02 

-1.08 

-1.28 

4 6 

-.18 

-07 

.06 

2.48 

2.64 

62 

4.45 

4.36 

-2.28 

5-53 

5-45 

a 

3-84 

3.52 

-.56 

4.12 

3.88 

66 

-5.25 

-5.65 

1.44 

-5-93 

-6.45 



Vertical displacements 

, in. 


22 

0.0288 

0.0283 

0.0282 

0.0133 

0.0247 

24 

.0290 

.0286 

.0288 

.0262 

.0282 

26 

.0283 

.0285 

.0283 

.0592 

.0336 

42 

.0875 

.0871 

.0876 

.2180 

.1887 

44 

.0880 

.0872 

.0870 

.2105 

.1870 

46 

.0895 

.0885 

.0880 

.1912 

.1830 

62 

•1577 

.1670 

.1735 

.2595 

-2515 

64 

.1660 

.1701 

.1702 

.2660 

.2530 

66 

.1926 

.1770 

.1627 

•3005 

.2605 


Condition (l): Vertical loads applied to points A and A' of end ring (see fig. 12(a)). No 

cutouts. Moment of Inertia of rings, 5*33 in.’’. 

^Condition (2) : Same as condition (l) except moment of inertia of rings increased by a 

factor of 5- 

°Condition (3): Same as condition (l) except vertical loads applied to points B and B' 

(see fig. 12(a)). 

^Condition (4): Same as condition (l) except with symmetrical cutouts in center hay 

(see fig. 12(a)). 

e Condition (5): Same as condition (4) except moment of inertia of rings Increased by a 

factor of 5* 
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(b) u^. circuit. 

Figure 6.- Electrical analogy for stringers and skin of a circular 

noncylindrical shell. 
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Figure 




(b) Displacements. 

7 .- Segment of a ring showing applied loads, internal forces 

and displacements. 





(b) Analog circuit for quadrant of ring. 

Figure 9*- Electrical analogy for a circular ring subjected to concentrated 
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(a) Flat rigid bulkhead. 


u t = e 



(b) Analogous electric circuit. 

Figure 11.- Electrical analogy for a flat rigid bulkhead. 
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(d) ut circuit. 
Figure 12. --Concluded 
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